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The cooling phase of Type I X-ray bursts observed with 
RXTE in 4U 1820-30 does not follow the canonical F oc T 4 
relation 
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ABSTRACT 

We analysed the complete set of bursts from the neutron-star low-mass X-ray binary 
4U 1820-30 detected with the Rossi X-ray Timing Explorer (RXTE). We found that 
all arc photospheric radius expansion bursts, and have similar duration, peak flux 
and fluence. From the analysis of time-resolved spectra during the cooling phase of 
the bursts, we found that the relation between the bolometric flux and the temper- 
ature is very different from the canonical F oc T 4 relation that is expected if the 
apparent emitting area on the surface of the neutron star remains constant. The flux- 
temperature relation can be fitted using a broken power law, with indices ^i=2.0±0.3 
and ^2=5.72±0.06. The departure from the F oc T 4 relation during the cooling phase 
of the X-ray bursts in 4U 1820-30 could be due to changes in the emitting area of 
the neutron star while the atmosphere cools-down, variations in the colour-correction 
factor due to chemical evolution, or the presence of a source of heat, e.g. residual 
hydrogen nuclear burning, playing an important role when the burst emission ceases. 
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1 INTRODUCTION 



Accreting neutron stars (NSs) in low-mass X-ray bina- 
ries (LMXBs) show episodes of unstable thermonuclear 
burning of hydrogen and helium on their surfaces (e.g. 
Lewin. van Paradiis. fc Taamlll993l ; IStrohmaver fc Bildstenl 



the expected surface area of a standard NS (i.e. ~400 km 2 ;, 
if a symmetric cooling phase is assumed jFrvxell fc Woosleyl 



19821: iBildstenll 19951: [Spitkovskv. Levin, fc Ushomirskv||200: 
Strohmaver fc Bildstenll2006l ) . This emitting area can be es- 



2; 



2006; iGallowav et alj|2008h . These episodes are called Type 



timated from the time-resolved spectral analysis during the 
cooling phase of the bursts. 



I X-ray bursts. Although the vast majority of the energy re- 
leased in LMXBs originates in the accretion process, capable 
of radiating ~300 MeV per accreted baryon, the unstable fu- 
sion process (<5 MeV per nucleon) during an X-ray burst 
can increase the luminosity of these objects by a factor of 
~10 in about 0.5-5 s. After that fast rise, the X-ray flux 
decreases again exponentially within 10-100 s. The typical 
total energy emitted during an X-ray burst is ~10 39 erg. 
Some of these bursts are strong enough to reach the Edding- 
ton luminosity. In those cases, the radiation pressure is high 
enough to trigger the expansion of the outer layers of the NS 
atmosphere, in a so-called Photos pheric Radius Expansion 
(PRE) burst |Basinska et al.lll984T ). The hypothesis that the 
bursts originate at the NS surface is supported by the fact 
that the inferred emission area during the bursts matches 
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The X-ray bursting LMXB 4U 
tracompact X-ray binary, with an 



1820-30 is an ul- 
orbital period of 



11.4 min jStella. White, fc Priedhorskvl 1 19871 ). consisting 
of a NS and a hydrogen-exha usted white-dwarf of 0.06- 
0.08 M m |Rappaport et al.|[l987l ). 4U 1820-30 belongs to the 
globular cluster NGC 6624, which is located at a distance 
of 7.6±0.4 kpc (|Kuulkers et al.ll2003l ). Th ermonuclear X- 
ray bu rsts in this LMXB were discovered bv lGrindlav et ail 
(1976). T he source also displays kilohertz quasi-periodic os - 
cillations (|Smale. Zhang, fc Whitell 19971 : IZhang et al.lll998l ). 
Until 2007 June, five Type I X-ray bursts -all PRE- 
were detected by the Rossi X-ray Timing Explorer, RXTE 
(|Gallowav et al.fl2008l). From time-resol ved spectra of three 
of these five bursts, Giiver et al.l (2010) found that, during 
the cooling phase, the bolometric flux was proportional to 
the temperature to the fourth power, and from this they 
estimated the mass and ra dius of the NS in 4U 1820-30. 
IStrohmaver fc Brownl (2002) also discovered a superburst in 
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this source, presumably due to the burning of a deep carbon 
layer. I n order to expla in the characteristics of these X-ray 
bursts. ICumminel (|2003u developed models of a NS accreting 
He-rich material. 

In this paper we study the complete sample of X-ray 
bursts, 16 in total, in 4U 1820-30 detected with RXTE. We 
explore the properties of these bursts and we focus on their 
behaviour during the cooli ng phase. Finally, we compare 
our re sults with those from IZhang. Mendez. fc Altamirand 
(|20 1 lT ) on the cooling phase of 4U 1636 -56, and with re- 
cent N S atmosphere models developed bv lSuleimanov et al.l 
<|201ll ). 



2 OBSERVATIONS AND DATA ANALYSIS 

We analysed all archival data from the RXTE Propor- 
tional Counter Array (PCA) of 4U 1820-30. These obser- 
vations cover the period from May 1996 to April 2011. 
The PCA detector consists of five collimated proportional 
counter units (PCUs) sensitive in the 2-60 keV energy 
range. In order to detect the X-ray bursts in the dataset, 
we generated 1-s time-resolution light-curves from the 
Standard- 1 mode (data packed in 0.125-s time-resolution 
bins without energy resolution) and we analysed each light- 
curve to search for three consecutive points that were 
8-cr above the average count rate of the observation. We 
have chosen this value to distinguish real bursts from fluctu- 
ations in the persistent emission. We confirmed each burst 
detection by inspecting all the Standard- 1 light-curves in the 
whole dataset by eye. We used Standard-2 mode data (16-s 
time-resolution and 129 energy channels) to produce X-ray 
colours of the sour ce, following the procedure described in 
IZhang et all |2009). We selected four energy bands: A — 
2.0-3.5 keV, B = 3.5-6 keV, C = 6.0-9.7 keV and D = 
9.7-16.0 keV interpolating in channel space, and we defined 
soft and hard colours as count rate ratios SC = B/A and 
HC = D/C, respectively. We also defined the X-ray inten- 
sity, /x, as the 2.0-16.0 keV count rate. The count rates in 
each energy band were normalized using the PCA data of 
the Crab source obtained close in time to our observations. 

To study in detail each detected X-ray burst, we used 
the PCA Event mode data, E_125us_64M_0_ls, where data is 
binned in 64 energy channels in the 2-60 keV energy range, 
with 125-/iS time-resolution. For the time-resolved analysis 
of the bursts we extracted spectra using the 64 channels in 
time-intervals going from, respectively, 0.125 s at the peak to 
1-2 s at the tail of the burst. We corrected each spectrum for 
dead-time using the methods supplied by the RXTE tearrQ 
and we generated response matrix files for the PCUs ac- 
tive during each observation containin g bursts. We f itted 
the data using XSPEC version 12.7.0 (|Arnaudlll996h . as- 
suming 0.5 per cent systematic error, as suggested by the 
RXTE calibration tearr0. For all time-resolved spectra (250 
in total) we fitted the data from 2.5 to 20 keV. An arti- 
fact in the binning scheme of the Event mo de data leads 
to zero counts in one channel per PCU unit |jahoda et al.l 



1 http://heasarc.gsfc.nasa.gov/docs/xte/whatsnew/deadtimc.html 

2 http: / /www. universe. nasa.gov/xrays/programs/rxte/pca/ 
doc/rmf/pcarmf-11.7/ 
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Figure 1. Colour-colour diagram of all RXTE observations of 
4U 1820-30. Grey points correspond to the whole set of obser- 
vations with 256-s time resolution. The colour of the persistent 
emission previous to the 16 X-ray bursts is indicated with red 
crosses. The colour of the persistent emission previous to the su- 
perburst is indicated with a blue star. Soft and hard colours are 
in Crab units. 
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Figure 2. Hardness-intensity diagram of all RXTE observations 
of 4U 1820-30. The symbols are the same as in Figure[T] Intensity 
and hard colour are in Crab units. 

1996). This is normally accounted for in the response ma- 
trix files. We found, however, that if only one PCU was on 
during an observation, that single channel with zero counts 
disrupted the fits. In our case, this happened in a hand- 
ful of spectra in which only PCU2 was on, and channel 
11 had zero counts. For this reason we ignored channel 11 
(E ~ 4 keV) in all our spectra. We ended up with 25 chan- 
nels in the 2.5-20 keV range for each spectrum. We extracted 
~100 s spectrum from the persistent emission before and 
after each burst and then used this, independently, as back- 
ground for the burst spectra, following the standard proce - 
dure in X-ray bursts analysis (e.g. iKuulkers et al.l [2002V 
We fitted the spectra using a blackbody model (bbodyrad) 
times a component that accounts for the interstellar absorp- 
tion towards the source (tbabs) fixing the hydrogen column 
density, Nn. The model for the interstellar absorption uses 
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Table 1. X-ray bursts in 4U 1820-30 detected with RXTE. 



Observation ID Starting time Duration SC * HC * ix * 

[MJD] [s] [Crab] [Crab] [Crab] 



+ 
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20075- 


01 


05 


00 


50570 731 1 6038 




oUUO f- 




n i 
U 1 




Atn n/ifii f^no 
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40017- 


01 


24 


00 


52794.73816848 
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70030- 


03 


04 


01 


52802.07560482 
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70030- 


03 


05 


01 


52805.89569741 


t 


90027- 


01 


0.3 


05 


53277.43859093 




94090- 


01 


01 


02 


54948.82125875 




94090- 


01 


01 


05 


54950.70282704 




94090- 


01 


02 


03 


54956.77474255 




94090- 


01 


02 


02 


54958.73999718 




94090- 


01 


01 


00 


54978.32150760 




94090- 


01 


04 


01 


54978.49491038 




94090- 


01 


05 


00 


54981.18730621 




94090- 


02 


01 


00 


54994.53419278 




94090- 


02 


01 


00 


54994.61302380 




96090- 


01 


01 


00G 


55624.88068005 




96090- 


01 


01 


020 


55626.77361987 



~26 


1.11(1) 


0.91(1) 


0.1263(3) 


>3000 ** 


1.08(4) 


0.80(4) 


0.112(1) 


-28 


1.19(2) 


1.00(1) 


0.0894(4) 


~26 


1.15(2) 


0.98(1) 


0.0914(4) 


~25 


1.19(1) 


1.01(1) 


0.1238(4) 


~23 


1.09(1) 


0.89(1) 


0.1317(5) 


~28 


1.15(1) 


0.98(1) 


0.1130(4) 


~30 


1.22(2) 


1.01(1) 


0.1012(4) 


~25 


1.26(2) 


1.09(2) 


0.0949(4) 


~31 


1.27(2) 


1.10(1) 


0.0912(4) 


~30 


1.33(2) 


1.11(2) 


0.1052(5) 


~29 


1.30(2) 


1.14(1) 


0.1048(5) 


~26 


1.29(2) 


1.12(1) 


0.1130(5) 


~21 


1.18(1) 


1.02(1) 


0.1855(6) 


~22 


1.19(1) 


0.99(1) 


0.1905(7) 


~25 


1.22(2) 


1.01(1) 


0.0940(4) 


~21 


1.11(1) 


0.83(1) 


0.1707(6) 



* Soft /Hard colours and intensity of the persistent emission just before the X-ray burst, in 

Crab units. 

** Super burst, see IStrohmaver fc Browij i2002| j. 
f X-ray bursts analysed bv lGuver*'et*ld f l|2010l ). 



the abundances o f I Wilms. Allen, fc McCravl (|2000h and the 
cross sections of IVerner et al. I (| 19961 ). The spectral model 
consists of two free parameters: the blackbody temperature, 
fcTbb, in keV units (where k is the Boltzmann constant) 
and the blackbody normalization, iVbb = (Rbb/dio) 2 , where 
-Rbb is the blackbody radius in km and dio is the distance 
to the source in units of 10 kpc. These parameters allow 
us to estimate the evolution of th e bolometric flux dur - 
ing the bursts using, e.g., eq. (3) of lGallowav et all (2008), 
F = 0.001076iVbb (fcTbb) 4 10" 8 erg cm" 2 s -2 . 

This procedure fails if the blackbody observed during 
the burst comes from the same blackbody that produces 
the persistent emission before the burst, because the dif- 
ference between two blackbo dy spectra is not a blackbody 
(van Paradiis fc Lewinlll986l ). However, in practice, the de- 
viations from a Planckian distribution become important 
only if the net emission is weak. Therefore, this issue may 
be important only at the very beginning or at the tail of the 
bursts. 

To test whether using the parameters of the fitting of 
the net spectrum instead of the total spectrum -which in- 
cludes modelling the persistent emission during the burst 
and subtracting the instrumental background- could affect 
our final results, we proceeded as follows: we fitted the per- 
sistent spectrum before each burst using an absorbed black- 
body plus a power-law model. Then, for the bursting time, 
we fixed the power-law index and normalization to the best- 
fitting values obtained for the persistent emission and we 
fitted only the blackbody parameters. The results from this 
procedure are consistent with those we obtained by fitting 
the net spectrum within errors, either using the persistent 
emission before or after each burst as background. To study 
whether possible variations in the persistent emission during 
the bursting time could change our results, we also fitted the 
spectra leaving free the power-law component. The results 
of this analysis are fully consistent within the errors with 



those found by fitting the net spectra. We also evaluated 
the contribution of the persistent emission to the net spec- 
tra during bursting time; we found that, for every spectra, 
the blackbody contributes at least 25 per cent of the total 
emission in the 2.5-20 keV energy band and in 75 per cent 
of the cases the blackbody contribution in this energy range 
is larger than 50 per cent. In the rest of this work we show 
the results from the standard procedure. 

In all spectra the number of counts per channel in the 
source and in the ba ckground is hi gh enough to justify the 
use of minimum x (|Pearsonlll900l ) to find the best-f itting 
parameters. Therefore, contrary to iGiiver et al.l |2010l ). we 
did not rebin the spectra in energy (we note also that the 
original binning scheme of the event-mode data that we used 
in the spectral analysis oversamples the intrinsic spectral 
resolution of the PCA by a factor of 1 to 3, depending on 
energy, and hence rebinning is neither required nor desir- 
able). For the sake of comparison, we also fi t ted al l the spec- 
tra using the C statistics defined by ICashl (|l979l ). which is 
the XSPEC implementation of the maximum-likelihood es- 
timate for the case of data following Poisson statistics. We 
found consistent results with both methods, and therefore 
here we report the results from the minimum \ 2 method. In 
the plots we show the 1-cr errors (A^ 2 = I). 



3 RESULTS 

3.1 X-ray burst detection and colour analysis 

We identified 16 short bursts (20-30 s), and the pres- 
ence of one superbur s t (>3 000 s), previously reported by 
IStrohmaver fc Brownl (|2002| ) in all the RXTE archival data. 
In Table [1] we summarize the main properties of the X-ray 
bursts. The starting time of the bursts was defined when a 
continuous increase in the PCA count rate persisted for at 
least 0.25 s within each detected burst (see Section 2), while 
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OBSID 96090-01 -01 -00G 




5 10 15 

Time [s] since MJD 55624.88068005 

Figure 3. The X-ray burst in 4U 1820-30 detected in ObsID 
96090-01-01-00G. The panels from top to bottom correspond to 
the Standard-1 mode (0.125-s time-resolution) X-ray light-curve 
(starting 1 s before the burst detection), the blackbody tempera- 
ture, fcTjjb, normalization, Nbb, bolometric flux, F, and reduced 
\ 2 for 23 degrees of freedom of the best-fit to the time-resolved 
spectra (starting from the onset of the burst), with 1-<t error-bars. 
The vertical dotted line indicates the touchdown point, which de- 
fines the starting of the cooling phase of the PRE burst. 



the ending time was denned as the time when the burst 
count rate returned to within 3-cr of the flux of the persis- 
tent emission before the onset of the burst. We also show 
in the table the burst duration -the difference between end- 
ing and starting times- and the soft and hard colours and 
the intensity of the persistent emission before each burst, 
derived from the Standard-2 mode data of the PCA. 

Figures [T] and [5J respectively, show the position of the 
source in the CD and HID before each burst. No bursts were 
detected for a persistent HC < 0.8 or for 1% > 0.2 Crab, 
where the s ource is bright and soft and th e accretion rate is 
likely high (|Hasinger fe van der Kli3ll989T >. 



before the onset of the burst. When the burst ignites, the 
light-curve shows a very high and narrow peak lasting for 
1-2 s. After that feature, the count rate starts to increase 
again reaching a maximum of ~5000 counts s _1 in 2-3 s, and 
finally it decays exponentially until the persistent count rate 
is recovered in 15-25 s. In the next four panels we present the 
evolution of the best-fitting parameters fclbb and iVbb, the 
bolometric flux derived from the blackbody model and the 
reduced x 2 residuals of the best-fitting model for 23 degrees 
of freedom, respectively. The distribution of the best-fitting 
X 2 of the 250 spectra is consistent with the expected x 2 dis- 
tribution. A Kolmogorov-Smirnov test gives a probability of 
0.09, which is not enough to disprove the hypothesis that 
the best-fitting x 2 are drawn from the expected x 2 distribu- 
tion for 23 degrees of freedom. During the spectral fitting, 
we considered tw o different values for the hydro gen column 
densit y, following ICostantini et al] (I2012T ) and ICiiver et al.l 
(|201Qh : iV H = 0.16 and 0.25 x 10 22 cm" 2 , respectively. The 
best-fitting spectral parameters obtained for these column 
densities are fully consistent with each other, because in 
the 2.5-20 keV energy band, absorption is not very sig- 
nificant. From now on, we show results obtained assuming 
iV H = 0.25 x 10 22 cm" 2 . 

At the beginning of the burst the blackbody temper- 
ature, feTbb, starts at ~3 keV and decreases rapidly, while 
the blackbody normalization, Nbb, which is proportional to 
the inferred emitting area, increases very fast. Then, fclbb 
starts to increase again as Abb decreases, until the touch- 
down point is reached. At this point, the A?bb is minimum 
and the temperature is maximum. Thus, we define this time 
as the start of the cooling phase. This common behaviour 
seen in all bursts indicates that they are all PRE. The high 
and narrow peak obser ved in the S tandard-1 light curves at 
the start of the burst ((Keck 2012]) is due to the fact that 
the temperature decreases enough to move the Planckian 
flux distribution to the soft band, outside the energy range 
where the PCA is sensitive dTawara. Havakawa. fc Kii|[l984l ; 
lin't Zand fc Weinberd[2f?Iol ). 

In order to study the properties of the whole sample 
of bursts, we analysed the behaviour of four parameters in- 
ferred from the X-ray analysis with respect to the position 
of the bursts in the CD. These parameters are: the temper- 
ature, fcTpcak, and flux, -Fpeak, at the touchdown point, the 
fluence, B, defined as the total energy per unit area emitted 
during the burst, B — J F(t)dt, and the decay time scale, r, 
obtained from the fit of an exponential function to the flux 
during the cooling phase of the burst. We divided the bursts 
in two groups, based on the hard colour of the persistent 
emission at the time of the burst: HC > 1.05 and HC < 1.05. 
The statistical analysis of both subsamples shows that both 
populati ons are indistinguisha ble through these parameters 
(see also lin't Zand et alj|20l3 ). 



3.2 Time-resolved spectral analysis of the X-ray 
bursts 

In Figure [3] we show the results of the fits to one of the X- 
ray bursts detected in 4U 1820-30. Because all the bursts 
in our dataset are very similar, we only show one of them 
as an example to illustrate the observed burst characteris- 
tics and the analysis process. In the top panel we present 
the Standard-1 high time-resolution light-curve, starting 1 s 



3.3 Flux-temperature relation during the cooling 
phase 

In Figure 2] we explore the relation between the bolometric 
flux, F, and the blackbody temperature, kTbb, during the 
cooling phase of all the busts. The diagonal lines in the plot 
represent lines of constant radius, R = 4, 6 and 8 km, re- 
spectively, assuming a distance d = 7.6 kpc to the source. 
Besides the fact that the dispersion is high at low flux levels, 
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Figure 4. Merged plot of the cooling phase of the 16 X-ray bursts 
of 4U 1820-30 detected with RXTE. The three diagonal lines 
are lines of constant radius, corresponding, respectively, to three 
blackbodies of _R = 4, 6 and 8 km, assuming a distance of 7.6 kpc 
to the source. We show the typical error bars for three different 
bolometric flux levels at the right side of the plot. 
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Figure 5. Same data as in Figure|4] rebinned into 25 points. The 
data are fitted by a broken power-law with indices v\ = 2.0 ± 0.3 
and V2 = 5.72 ± 0.06 above and below the break, respectively. 
The break is at kT hh = 2.43 ± 0.02 keV. 



all bursts follow the same trend in this diagram (the spread 
in the plot is consistent with the statistical errors of the flux 
and temperature). Since the flux-temperature relation of all 
bursts is consistent with being the same, we binned the data 
sorted by flux into 25 points, and we calculated the weighted 
average of the flux and temperature within each bin. We also 
binned the data sorted by temperature, and we used differ- 
ent binning factors, but in all cases we obtained consistent 
results. In Figure we plot the flux-temperature relation 
for the rebinned data. From this plot it is apparent that 
the cooling phase of the X-ray bursts in 4U 1820-30 does 
not follow the expected F ex (fcTbt,) 4 relation for an ideal 
blackbody of constant emitting area. Fitting the data with 
a single power law, F oc (kTbb)", gives ^=4.95±0.03, but 
the fit is not good, with x 2 =499 for 23 degrees of freedom. 
A broken power-law fits the data better (although formally 
it is not a good fit), with x 2 — 44 for 21 degrees of freedom, 



resulting in power-law indices before and after the break 
z/i=2.0±0.3 and i/2=5.72±0.06, respectively, and a break at 
feT bb =2.43±0.02 keV (see Figure |5J|. 

We also studied the flux-temperature relation dividing 
the bursts into two groups for HC > 1.05 and HC < 1.05, as 
in the previous subsection. After rebinning both sets of data 
into 25 points, we plotted their respective cooling phases 
and confirmed that both curves are consistent within errors, 
which lends support to the idea that the full sample of bursts 
in 4U 1820-30 is homogeneous. 



4 DISCUSSION 

We present, for the first time, a homogeneous analysis of the 
cooling phase of all the thermonuclear X-ray bursts detected 
with RXTE in the ultracompact LMXB 4U 1820-30. The 
sample contains 16 PRE bursts with a duration of 20—30 s, 
all of them detected when the source was in the hard state 
or in the transition from the hard to the soft state, at rel- 
atively low luminosity. The touchdown temperatures, bolo- 
metric peak fluxes, fluences and time-scales of the exponen- 
tial cooling-tails are consistent with being the same in all 
bursts in the sample, supporting the idea that the bursts 
in 4U 1820-30 are very similar between them. Furthermore, 
contrary to what is expected for a constant-radius black- 
body, the cooling phase of the X-ray bursts do not follow 
the F oc (fcTbb) 4 relation. In fact, the relation between the 
bolometric flux and the blackbody temperature is fitted bet- 
ter by a broken power-law of indices 2.0T0.3 and 5.72T0.06 
than with a single power law. The observed deviation from 
the F oc (fcTbb) 4 relation complicates the use of the cooling 
phase of X-ray bursts to infer the mass and radius of the NS 
in 4U 1820-30. 

The global trend that we found for the cooling phase 
of the X-ray bursts in 4U 1820-30 is similar to that of the 
hard non-PRE bursts in 4U 1636-53 (Zhang et al. 2011). In 
both cases, the bolometric flux decays faster than what is 
expected for a blackbody with constant radius. This is re- 
flected by the value of the power-law index (v2 = 5.72T0.06) 
that best-fits the cooling phase of the bursts below kTbb ~ 
2.43 keV (see Figure[SJ), and also in the best-fitting value for 
a single power law (y = 4.95 ± 0.03) to the whole cooling- 
phase data. Taking into account the whole set of bursts in 
4U 1820-30 observed with RXTE, there is no region of the 
F vs. fcTbb diagram that can be fitted by a power-law with 
an index equal t o 4. O ur results are at variance with those 
of [G iiver et all (|2010l ). who analysed five bursts from 4U 
1820-30 observ ed with RXTE (ind icated with a f symbol 
m our Table Q) . iGiiver et al.l (|201Ch showed that during the 
cooling phase, three of these bursts followed the relation 
F oc (fcTbb) 4 , whereas they discarded the other two bursts 
because, according to their own explanation, these bursts 
showed significant and continuous fluctu ations in the appar- 
ent radius. The three bursts analysed bv lGiive~ et all (|2010h 
are included in our Figure [4] where it is apparent that none 
of the bursts in 4U 1820-30 follows the F oc (fcT bb ) 4 relation, 
even in the narro w flux range (1.0-6. Ox 10 -8 erg cm~ 2 s _1 ) 
that lGiiver et all used. We fitted a power-law to these three 
bursts in this flux range obtaining an index of 4.8±0.16 (1-a 
error; \ 2 = 15 for 16 degrees of freedom), which is signif- 
icantly different from 4. Fitting a power-law to the whole 
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sample of bursts in the same flux range gives an index of 
4.68±0.08 (1-<T error; \ 2 = 162 for 98 degrees of freedom). 
We can not explain the dif ference between our results and 
those of lGiiver et"al] I^PIOD . given that, as described in Sec- 
tion 2, we analysed the data as they did. We note, however, 
that from our analysis, all bursts in 4U 1820-30 observed 
with RXTE are consistent with following the same trend in 
the flux-temperature diagram during the cooling phase, in- 
cluding the three bursts studied by iGiiver et all but also 
those two bursts that they discarded. 

The departure from the F oc (fcTbt,) 4 relation during the 
cooling phase of the X-ray bursts in 4U 1820-30 could be due 
to three different -but not mutually exclusive- reasons: (i) 
the non-constancy of the true emitting area of the NS dur- 
ing this phase; (ii) changes of the colour-correction factor, 
f c — Tbb/Teff, which accounts for hardening of the spectrum 
due to (Compton) electron scatte ring in the NS atmosphere 
ijLondon. Taam. fc Howardll 198(f ) during the cooling phase; 
(iii) the presence of a source of heat that increases the ob- 
served temperature, e.g. residual hydrogen nuclear burning 
during the cooling phase of the bursts. A final possibility is 
that the real underlying spectra are not blackbodies. 

If f c remains constant during the cooling phase, varia- 
tions in R oc . could be explained by changes in the true emit- 
ting area of the NS. In that case, if a hotspot is formed away 
from the rotation p ole, it should be possible to ob serve burst 
oscillations (see e.g. lStrohmaver fe~ Bildstcn 2006). However, 
no oscillations have been detected in t he bursts of 4U 1820- 
30 (|Gallowav et al.ll2008l ; IWattdl2012h . Alternatively, if the 
emitting area does not change, variations in f c could ex- 
plain these differences. In terms of the spectral parameters, 
the colour-correction factor can be expressed as: 



fc 



Tcff 



Roo 
Rbb 



\ 



Ro 



F 
trT' 1 , 



R(l 



(1) 



where Roo is the NS radius seen by a distant observer, z is 
the gravitational redshift at the surface of the NS, d is the 
distance to the source, R the radius of the NS, and a the 
Stefan-Boltzmann constant. 

In Figure[6]we show the plot of the inferred f c as a func- 
tion of the luminosity in Eddington units. In the plot, we 
obtained the colour-correction factor from the spectral pa- 
rameters, assuming R — 14 km, z — 0.19 and d = 7.6 kpc. To 
estimate the luminosity in Eddington units we assumed that 
the flux reaches the Eddington limit at the touchdown point 
in each burst, and further that the emission is isotropic. 
Hence, we averaged the peak flux of the 16 bursts, obtaining 
-Ffedd = 6 x 10~ 8 erg cm -2 s - 1 . In the plot we also sho w the 
theoretical models of f c from ISuleimanov et al.l l|201l[) that 
correspond to atmospheres with three different chemical 
compositions: pure hydrogen, solar abundance and pure he- 
lium, respectively. Starting at near-Eddington luminosities, 
as the flux decreases, the inferred f c slightly decreases until 
it reaches a minimum value / c ~1.4 close to L/L^dd = 0.6. 
After that, the inferred colour-correction factor starts to in- 
crease again, doing so even faster for L/L-^dd < 0.2, in the 
cooling tail of the bursts, close to the level of the persistent 
emission, where the errors in the fitted spectral parameters 
become important and the models are less accurate. Besides 
the fact that the models do not fit the data well, which pre- 
vents us from constraining the atmosphere composition, our 




Figure 6. Average colour-correction factor. / c , as a function 
of luminosity in Eddington units for the 16 X-ray bursts in 
4U 1820-30. The points correspond to the whole sample of bursts 
in 25 bins, using the data in Figure [5] and eq. JTJ. Dashed lines 
correspond to models of f c for three d ifferent chemical compo - 
sitions of the neutron-star atmosphere ( Sulcimanov et al.ll20l JV 
indicated in the legend. 



results are incompatible with a solar-type atmosphere be- 
cause the data do not show any dip at L/L^dd < 0.2. At the 
beginning of the cooling, the data are consistent with mod- 
els with He-rich composition, which agree s with the chem- 
ical c omposition of the accreted material l|Rappaport et all 
1987) deduced from the properties of the white-dwarf com- 
panion. Moreover, the shape of / c in this plot is similar to 
that of the hard non-PRE bursts in 4U 1636-53 (Zhang et 
al. 2011). 

We note also that the distribution of fcTbb at constant 
bolometric flux levels broadens close to the tail of the bursts 
(Figure |4]). This could be due to differences in the under- 
lying emission coming from the NS surface. This emission 
can be due to persistent accretion playing an important role 
when the burst emission ceases, or to the burning of a small 
residual fraction of hydrogen, not fully burned before touch- 
down, during the tail of the bursts. Hydrogen burning is 
slower than helium burning because the /3-decay within the 
nuclear chain is driven by the weak force, thus having a 
longer time-scale than the triple-a process. The possibility 
that this source o f heat operates d uring the cooling phase of 
the X-ray bursts jCumming 2003), could explain the power- 
law index larger than four that best-fits the relation between 
the flux and the blackbody temperature. 
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